function Uout = hamil2Q(Uin,np,fieldeps,E0,muop)

%        Hamiltonian operatioon 

%        phi = H psi

        Uout = 0;
        Uout = (fieldeps*muop+diag(E0))*Uin;
        end